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ABSTRACT 

Context. Polarimetry is one of the keys to enhanced direct imaging of exoplanets. Not only does it deliver a differential observable pro¬ 
viding extra contrast, but when coupled with spectroscopy, it also reveals valuable information on the exoplanetary atmospheric com¬ 
position. Nevertheless, angular separation and contrast ratio to the host-star make for extremely challenging observation. Producing 
detailed predictions for exactly how the expected signals should appear is of critical importance for the designs and observational 
strategies of tomorrow’s telescopes. 

Aims. We aim at accurately determining the magnitudes and evolution of the main observational signatures for imaging an exoplanet: 
separation, contrast ratio to the host-star and polarization as a function of the orbital geometry and the reflectance parameters of the 
exoplanet. 

Methods. These parameters were used to construct polarized-reflectance model based on the input of orbital parameters and two 
albedo values. The model is able to calculate a variety of observational predictions for exoplanets at any orbital time. 

Results. The inter-dependency of the three main observational criteria - angular separation, contrast ratio, polarization - result in a 
complex time-evolution of the system. They greatly affect the viability of planet observation by direct imaging. We introduce a new 
generic display of the main observational criteria, which enables an observer to determine whether an exoplanet is within detection 
limits: the Separation-POlarization-Contrast diagrams (SPOC). 

Conclusions. We explore the complex effect of orbital and albedo parameters on the visibility of an exoplanet. The code we developed 
is available for public use and collaborative improvement on the python package index, together with its documentation. It is another 
step towards a full comprehensive simulation tool for predicting and interpreting the results of future observational exoplanetary 
discovery campaigns. 
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1. Introduction 


Observational technologies for detecting the light reflected from 
an exoplanet are reaching a level of precision that makes direct 
imaging of exoplanets a realistic possibility in about the coming 
decade. Because of the very challenging contrast ratio between 
exoplanets and host-stars in the optical, there is great interest 
in differential methods, such as polarimetry, to deliver an ex¬ 
tended reach to imaging instruments (Hough et al. (|2006 1 , Keller 


( |2006| >, |Schmid et al.| ([20061). Built upon previous detections 
of polarized signals from exoplanets Berdyugi na et al.| ([2 008 , 
Wiktorowicz (|2009 1 , Lucas et al. (2009l, Berdyugina et al. 


( |201 1| >, the polarimetric imaging modules that are currently be¬ 
ing integrated on the 10m class telescopes (SPHERE-ZIMPOL 


et al. 

< 

et al. 

4 


2012), and VAMPIRES on SUBARU Telscope Norris 


2015)) promise to achieve a differential polarized con¬ 
trast down to approx 10~ 6 in the visible or near- infrared. It is 
therefore an opportune moment to produce detailed predictions 
for exactly how the expected signals should appear, which will 
be of critical importance for the designs and observational strate¬ 
gies of these instruments. We aim at accurately determining the 
magnitudes and evolution of the main observational signatures 
as a function of the basic parameters of the exoplanetary system: 
the star-exoplanet orbital parameters and the optical properties 


of the planet. We incorporate the complexity arising from effects 
such as polarization based on Rayleigh scattering. Additionally 
and as a first-order observability estimator, the absolute flux in 
photon per unit time and surface received from the exoplanet tar¬ 
get is computed assuming a black-bodied star. We finally provide 
a querier and parser of the http: //exoplanet. eu/ exoplanet 
database for searching and importing any star-planet target. 


Several models were already developed to predict the light 
signature of an exoplanet as it should appear to an observer, for 
example see |Cahoy et al.| ( |2010| ), |Buenzli & Schmid| ( |2009| l, or 
for the case of Earth-like planets with varying atmospheric pa¬ 
rameters, see Karalidi & Stam (2012]) and Zugger et al. ( 2010[ ). 
However, to the best of our knowledge, the inter-dependency of 
the main observational criteria is never taken carefully into ac¬ 
count. The relative evolution of angular separation, polarization, 
and contrast ratio is of critical importance because their respec¬ 
tive maxima do not occur at the same time in the general case: 
this makes the best-case scenario for detection very unlikely. We 
correct for this by showing the relative evolution of the angu¬ 
lar separation, contrast ratio, and polarization as a function of 
the orbital and reflectance parameters of the planet. The findings 
resulting from this integrated treatment highlight dependencies 
that are much more complex than previously reported. The three 
main observational criteria - angular separation, contrast ratio, 
and polarization - in general do not exhibit a strong positive cor- 
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relation and must be analyzed separately to determine the direct 
visibility of a potential planet target. 

2. Model 

We developed a polarized reflectance model. Its algorithm relies 
on the orbital parameters of the exoplanet and the linear combi¬ 
nation of two albedo values; it is briefly described below. 

As our prime observables are in any event differential be¬ 
tween the star and exoplanet (such as the contrast ratio) or in¬ 
herent to the planet (such as the polarization fraction), absolute 
fluxes are not described in detail. A simple black-body model¬ 
ing of the host stars is used to give a first-order idea on the ab¬ 
solute fluxes that are expected from a exoplanet. For the sake of 
conciseness, “polarization” is used instead of “polarization frac¬ 
tion”. 

2.1. Orbital parameters 

Each exoplanet orbit is defined by the usual parameters listed in 
Table [T] Note that the shape of the orbit as seen by the observer 
is solely defined by e, i, and u> p . 



Fig. 1. The geometric albedo Ag^ ay / cig /,arising from single¬ 
scattering as a function of the Rayleigh single-scattering albedo 
w R . 


Symbol 

Name 

Unit 

T 

Period 

days 

tp 

Time at periapsis 

Julian date 

a 

Semi-major axis 

AU 

e 

Eccentricity 

0 

i 

Inclination to the observer 

degrees 

OJp 

Argument of periapsis 

degrees 

n 

Longitude of ascending node 

degrees 


Table 1 . Orbital parameters for an exoplanet. /=90° corresponds 
to an edge-on orbit (transiting), ;^90° corresponds to a retrograde 
orbit. 


The phase angle of the exoplanet and the distance between 
the star and the planet are obtained from these orbital parame¬ 
ters. The phase angle ais defined as the vector angle between the 
star, an exoplanet, and the observer, as seen from the exoplanet. 
It is 0° when the planet is at full phase (superior conjunction) and 
180° when the exoplanet is at new phase (inferior conjunction, 
transiting). 

2.2. Rayleigh scattering polarization 

For the sake of conciseness, in the following discussions we use 
the quantity called phased albedo, which is defined as the prod¬ 
uct of the phase function and the geometric albedo of a planet, 
for a given wavelength: 

A a (a) = A g ■ (p(a) = (1) 

for an incident flux nh\ with 1(a) the emerging flux from a body 
at the phase angle a. It represents the fraction of the incident 
irradiance that is reflected by the planet when it is seen at phase 
angle crby the observer, so that E PlaneU ou , = A a (a) ■ E P i aneU in , 
with E being the spectral irradiance arriving at or leaving the 
planet. 

The polarization induced by an exoplanet on a reflected 
beam of light is described by many different models, such as 
Mie, Rayleigh, or Raman scattering. It is assumed here that sin¬ 
gle Rayleigh scattering is the predominant source of polarization 


in planetary atmospheres. This is a good approximation as long 
as Mie scattering (and especially the primary rainbow polariz¬ 
ing effect of liquid droplets in clouds, see Bailey ( |2007 1 remains 
lower than Rayleigh scattering. This scenario is usually achieved 
for wavelengths shorter than approximately 1 pm. which also cor¬ 
responds to the most favorable wavelengths for Rayleigh polar¬ 
ization measurements of scattered light, as it is function of A~ 4 . 
Integrated over the stellar disk, the flux from the star can be con¬ 
sidered to be unpolarized (Kemp et al. ( |1987[ ). Hence, only the 
light from the exoplanet carries polarization. 

The results of |Madhusudhan & Burrows] ( |2012j > for the 
Stokes parameters Q/I and (/// were adopted as input data to 
carry out more polarization and reflectance computations. Their 
model assumes an unresolved semi-infinite homogeneous atmo¬ 
sphere (hence cloud free), dominated by Rayleigh scattering. 
The atmosphere is assumed to be spherical, which is a satis¬ 
factory assumption when the planetary rotation remains slow. 
Given the strong depolarizing effect of multiple scattering, they 
used the single-scattering albedo as a unique reflectance param¬ 
eter from which a geometric albedo arising from the 

Rayleigh single-scattering albedo was calculated; their result is 
reproduced in Figure [I] 

Stokes-V is found to be zero for Rayleigh scattering in planet 
atmospheres. Hence, the degree of polarization PzR ay i e igh is de¬ 
fined by 


PZ-Rayleigh 




+ u l 

out ^ ^ out 


( 2 ) 


and it depends on the phase angle aof the planet and 
the single-scattering albedo wr , which represents the amount 
of absorption versus scattering in a given atmosphere. The 
existing literature provides solutions to several phase func- 
Madhusudhan & Burrows (2012)> provided an analyti- 


tions. 


cal solution for scattering models for Lambertian, Rayleigh, 
isotropic, and asymmetric scattering. Zugger et al. (2010) ad¬ 


dressed the liquid surface scattering model. Some phase func¬ 
tions are reproduced in Figure [2] Given the symmetry of the 
phase angle aalong the orbit, only a = [0,180]° values are repre¬ 
sented. The illuminated fraction values are also plotted for com¬ 
parison. 
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Fig. 2. Phase functions 0as a function of phase angle of or sev¬ 
eral scattering models. The illuminated fraction is also displayed 
(thick line). The Rayleigh-scattering model corresponds to a 
vectorial Rayleigh phase matrix. 


Rayleigh and Lambertian surface models follow the same 
general S-curve shape between zero flux and full flux. A 
Lambertian or Rayleigh planet is faint at small abecause the il¬ 
luminated fraction // is small, hence fewer photons are reflected. 

2.3. Polarized and unpolarized albedo 

It is now apparent that clouds and condensates are very com¬ 
mon in planet atmospheres. An atmosphere solely described by 
a single-scattering albedo could only take into account Rayleigh 
scattering of the atmosphere in the optically thick case, which 
is quite restrictive. Indeed, clouds, rocky surfaces or liquid sur¬ 
faces do not follow Rayleigh scattering, but they often strongly 
contribute to the reflectivity of the exoplanet. 

In this section we describe a novel approach to the re¬ 
flectance properties of an exoplanet atmosphere to include ef¬ 
fects of a greater atmospheric variety. 

We assumed that any non-Rayleigh single-scattering from 
the atmosphere follows Lambertian scattering. This assumption 
is discussed below. A planet is better defined by two combined 
albedo quantities: 

- a total geometric albedo A g j ota i\hdA represents the total frac¬ 
tion of incoming irradiance reflected by the exoplanet, taking 
into account all scattering processes, 

- a Rayleigh single-scattering albedo wr that represents the 

fraction of incoming irradiance that is scattered by the atmo¬ 
sphere according to Rayleigh single-scattering alone, from 
which arises a Rayleigh geometric albedo Only 

this part of reflected light carries a defined degree of polar¬ 
ization in this model. This albedo can be approximated by 
the single-scattering albedo of the exoplanet multiplied by 
the fraction of Rayleigh-scattered emerging light. 

These albedo values can be added linearly so that 

Ag'Total — Ag,Lambertian ^g,Rayleigh- ( 3 ) 


Following this, we can express the outgoing irradiances us¬ 
ing the phased albedo A a \ 

Act,Rayleigh ~ $Rayleigh ' A g,Rayleigh (4) 

A a,Total = 0 Lambertian ' ( A g,Total ~ Ag,Rayleigh) 

+ 0 Rayleigh ' Ag,Rayleigh ■ (5) 

The total emerging spectral irradiance of the exoplanet 
I Total,out including thermal emission is 

lTotal,out E'Rayleigh,out ELambertian,out Ej' hernial,out- ( 6 ) 

Note that this thermal emission is much weaker than the reflected 
light from the host-star unless the exoplanet is a hot Jupiter and 
the observing wavelength lies in the infrared. Below, we assume 
that E ThermaUmt is negligible for A obs < 1 pm. 

The different Stokes parameters and polarization degree can 
be easily obtained with 


QH 


I Rayleigh,o 


-)' E in -A, 


a,Rayleigh 


u/i = 


I Total ,out 
' Ei ’• ' Aa -Rayleigh 


It otal,t 


vi I = 0 

PZ-Total ~ 


PZ-Raylgeigh 


1 + 


(7) 

( 8 ) 

(9) 

( 10 ) 


where QlI Ray ieigh,out and U / I Ray ieigh,o,,t are the Stokes ratios from 
the Rayleigh-scattering polarization section of this model; they 
are used to calculate PzR ay i g eigh using Equation[2] 

According to previous sections, the following complete re¬ 
lation for the exoplanet reflected irradiance, given the host-star 
surface irradiance, is 


E 


Planet,Distance 


- Eg tar. 


Surface 


P* )2 _ f Planet ^ 
rorbit rl 


( 11 ) 


where R are radii. The observer-exoplanet distance is here ap¬ 
proximated by the distance observer-star d. The contrast ra¬ 
tio CRbetween a planet and its host-star is then obtained from 
Equation [IT] 

Cr = (fllyyyLf . A a . ( 12 ) 

rorbit 

We highlight here that the polarization of the exoplanet 
Pzrotai an d the contrast ratio CRdo not depend on the incom¬ 
ing irradiance from the host-star, they are intrinsic to the planet. 
However, they are strongly dependent on the observation wave¬ 
length /lobsthrongh Rayleigh single-scattering and geometric 
albedo, even though for clarity of notation, A was not explicitly 
written in the relations. 

Thanks to the separation of the contributions from the dif¬ 
ferent scattering processes on an exoplanet, an original model 
for planetary polarized reflectance has now been set up. It yields 
the phased albedo A a , which indicates the fraction of reflected 
irradiance from an observed planet, as a function of geometric 
albedo A g , phase angle a , and Rayleigh single-scattering albedo 
Wr. While crcan be easily determined as described in the orbital 
computations part of the model, the albedo parameters A,and 
Wfidepend on factors too numerous to be modeled here: they 
therefore need prior computation. If A ? and vi’ A >are obtained sep¬ 
arately, the albedo and phase functions previously defined are 
successful in describing the emerging irradiance that is reflected 
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by a planet. Its most accurate results are obtained for wave¬ 
lengths shorter than approximately 1 pm unless thermal emission 
and Mie scattering are also added to the model. 

We assumed a polarized-reflectance model that incorporates 
Rayleigh and Lambertian scattering as two linear scattering phe¬ 
nomena. The approximation made is that any scattering that 
does not relate to Rayleigh single-scattering is assumed to be 
Lambertian. As a consequence, a beam of light undergoing 
multiple-Rayleigh scattering or a combination of Lambertian 
and Rayleigh scattering is assumed to be multiple Lambertian 
scattered. This approximation mainly affects the angle at which 
a photon emerges from the exoplanet and in turn the phase- 
function of integrated emerging light. In our model, this inte¬ 
grated phase function is a linear combination of Rayleigh and 
Lambertian phase functions. Adding these second-order scatter¬ 
ing phenomena to the computation of the integrated phase func¬ 
tion adds terms to its computation: a multiple Rayleigh phase 
function term and a coupled Rayleigh-Lambertian phase func¬ 
tion term. We assumed that these second-order phase functions 
are similar to the Rayleigh and Lambertian phase functions, 
which differ only slightly (see Figure [2|. Furthermore, we as¬ 
sumed that the weight of these additional phase functions in 
the computation of the integrated phase function is smaller than 
Rayleigh single-scattering and Lambertian single or multiple 
scatterings. Note also that the lower the albedo values, the less 
likely multiple-scattering becomes and the smaller the errors in¬ 
duced by this approximation. 


The results were benched-marked against several other mod¬ 
els providing photopolarimetric curves as a function of one or 
several orbital parameters, such as models developed by Buenzli 
& Schmid ( |2009|l, |Madhusudhan & Burrows| ( |2012| l, Fluri & 
Berdyugina ( 2010[ >, and |Zugger et al.| ( 2010[ >. The locations of 


the polarization peaks and minima were reproduced with very 
good agreement; they mostly depend on the phase functions and 
the processing of orbital parameters. The intensity of the polar¬ 
ization peaks were reproduced with good accuracy; they mostly 
depended on the fine-tuning of the two albedo values. More 
specifically, the shift of the polarization peak to phase angles 
greater than 90° reported in the last bench-mark reference was 
also observed. 


2.4. Model completeness 


The calculation code that implements this model takes up to 11 
input parameters and N orbital positions for which the calcula¬ 
tions are performed. It generates N-element vectors for up to 11 
output parameters, listed in Table [2] Table [3] shows the mapping 
between the input and output parameters, where the input pa¬ 
rameters are listed as column headers and output parameters are 
row headers of Tabled 

In the way the code is written, the only mandatory input is 
the eccentricity, the inclination, and the argument at periapsis. 
When compared with exoplanetary detections reported in the lit¬ 
erature, these three orbital parameters are rarely known with any 
precision. This is currently the case for 132 exoplanets out of 
the 1790 («7.4%, source: http: //exoplanet. eu/). However, 
where more complete data exist, reasonably good estimates can 
be made. 


Output 

Name 

Unit 

a 

Phase angle 

degrees 

fr 

Period fraction since last periapsis 

% 

V 

True anomaly 

degrees 

r orbit 

Star-planet distance 

AU 

e 

Angular separation with star 

arcsec 

date 

Julian Date 

JD 

CR 

Contrast ratio with star 

0 

Pz 

Polarization degree 

% 


Polarization angle 

degrees 

(PNorth 

Angle with the north 

degrees 

Uorb 

Apparent orbital velocity 

mas/hour 


Table 2. Output parameters from the calculation code. Their pro¬ 
cessing depends on the availability of input parameters, refer to 
Table® ‘ ‘mas” stands for milli-arcsecond. 



e, i, (x) p 

a 

d 

T 

tp 

Rp, A g 

W R 

n 

a , f T , v 
r or bi, l ,0' 

X 








date 1 

X 



X 





date 1 

X 



X 

X 




T orbit > @ 

X 

X 

X 






~cr' 

X 





X 



“ CR 71 

X 

X 




X 



Pz, Oft 1 

X 





X 

X 



X 





X 

X 

X 

< PNorth 

X 







X 

Vorb 

X 

X 

X 

X 




X 


1 relative values 

2 absolute values 

Table 3. Mapping between the input (columns) and output 
(lines) parameters of the model. An “x” indicates that the given 
input parameter(s) is required for computing that output parame¬ 
ter)s). Example: calculating the relative angular separation (line 
4) and the absolute contrast ratio (line 7) only requires six input 
parameters: e, i, a> p , a, Rpi ane t , and A g . 


3. Results 


3. 1. SPOC diagrams 

An important tool introduced here is the Separation- 
POlarization-Contrast (SPOC diagram) which presents all 
the useful information to enable an observer to evaluate whether 
a planet target is observable at given instrumental limits and 
how this signal will evolve with orbit and time. Perhaps most 
importantly, in the event of a detection, SPOC diagramsprovide 
a powerful mechanism to constrain exoplanet properties given 
observational imaging data. 

Figure [3] shows an example of a SPOC diagramapplied 
to Alpha Centaury Bb, an Earth-mass planet discovered in 
(201 2[> that was later debated 


2012 

d20l3 


Dumusque et al. 


Its probable or 


Hatzes 

rital parameters are shown in Table [4] 


The exoplanet was discovered with the radial velocities method, 
hence its inclination and argument of ascending node are un¬ 
known. For i, values in [15, 35, 60, 85]°were explored, whereas 
Qwas neglected because it does not change the shape of the or¬ 
bit: it only defines a rotation of the orbit locus as seen by an 
observer along the line of sight (i.e., it is only useful for project¬ 
ing the exoplanet location around the star onto RA-DEC axes). 
Its radius was calculated from its mass -sin(i) = 0.0036 • Mj upiu , r 
assuming the same density as Earth, which leads to M = [4.42, 


4 



































































Guillaume Schworer & Peter G. Tuthill: Predicting exoplanet observability in contrast, separation and polarization 


2.00, 1.32, 1.15] M Earlh and R = [1.64, 1.26, 1.10, 1.05] R Earth 
for the four previous i. In this example, the light from objec- 
tAlpha Centaury A received by the exoplanet is neglected be¬ 
cause it is roughly 2.10 5 fainter than that from Alpha Centaury 
B. Reflectance parameters were chosen to be similar to those of 
a Venus-like planet (A g =0.61. w p =0.85) in V-band. 


Set 

T [day ] 

Up [°] 

a [AU] 

e 

t„ [ JD] 

Cen a $ 

3.2357 

0 

0.04 

0 

55280.17 

Cen a $M 

3.2357 

246 

0.04 

0.34 

55282.53 


Table 4. Two sets of probable orbital parameters for Alpha 
Centaury Bb, see section 9 of the supplementary information 
of Dumusque et al.j(j2012|. The inclination is choosen from [15, 
35, 60, 85]°, which leads to M = [4.42, 2.00, 1.32, 1.15] M Earth 
and R = [1.64, 1.26, 1.10, 1.05] R Ear rh respectively (assuming 
Earth density). 


In Figure [3] Cen a 0.34 is shown with i— 15°. This figure dis¬ 
plays the polarization encoded with color, while the two axes 
plot the separation against contrast ratio; finally, a set of differ¬ 
ent markers give temporal information. The first striking feature 
of the plot is that the exoplanet does not run a “back and forth” 
locus in the separation-contrast ratio phase diagram. This is be¬ 
cause the apparent symmetry of the orbit is broken. When the 
orbital parameters /and care non-zero and to p is different from 
[0,90,180,270]°, the geometry of the orbit as it appears to the 
observer has no symmetry. In other words, the star is not at a fo¬ 
cal point of the apparent elliptical locus of the orbit. Note that the 
semi-major axis odoes not play a role in the shape of the orbit; 
it acts as a simple scaling factor. The photon flux on the right- 
hand side y-axis is given in photons per collecting-area per hour, 
integrated over V-band. It was calculated assuming that Alpha 
Centaury B is a black-body, T e ff = 5214 and M v = 1.33. The 
collecting area is 51.7/n 2 , which corresponds to a diameter of 8.2 
meter with a central obstruction of 2% (same as VLT telescopes). 

Figure [3] shows that the lowest and highest contrast ratios 
are reached close to the maximum and minimum phase angle 
(black and white disks). This corresponds to the inferior and su¬ 
perior conjunction of the exoplanet with its star. We note that 
the periapsis and apoapsis are not reached at these phase angle 
extrema because of the non-null value of a> p , which re-orients 
the orbit with respect to its host star. Polarization reaches two 
maxima near a=90°(the limit between the shaded and the bright 
area in the bottom left subplot of the same figure). This is linked 
to the fact that polarization is observed from Rayleigh scatter¬ 
ing, which has a peak polarization for 90°. However, in the case 
of non-entirely Rayleigh planets, the peak polarization occurs at 
slightly different phase angles. 

An even more interesting diagram is shown in Figure [4j It 
represents the apparent orbital motion in mas per hour plotted 
against the angular separation as x-axis and the polarized con¬ 
trast ratio as y-axis. The apparent orbital motion is found to be 
an important factor to take into account here because the period 
of Alpha Centaury Bb is only 3.24 days: it reaches 1.8 to 3.5 mas 
per hour, which limits the observer to a few hours of exposures, 
depending on the plate scale of its detector. The photon flux is 
extremely high as a result of the brightness of the star (1.3mag 
in V) and the large bandwidth of the V filter (550 ± 44nm). The 
SPHERE-ZIMPOF performance curve for lh exposure time us¬ 
ing a V-filter with the double-difference polarization calibration 
was added to the plot (see Fig. 4 in Roelfsema et al. ( |2014| )). 


4.3e+08 


3.8e+08 

3.4e+08~ 

C 

o 

o 

3.0e+08 q. 

2.7e+08 

2.4e+08 

2.1e+08 


0.020 0.025 0.030 0.035 0.040 

Angular separation [arcsec] 

I I Degree of polarization f%1 | I 

24.0 24.4 24.8 25.2 25.6 26.0 26.4 26.8 

Fig. 3. SPOC diagramfor Ce« ff> o .34 with i= 1 5°, see Table [4] as¬ 
sumed to be Venus-like regarding its reflectance parameters. It 
represents the polarization values as gray scale, plotted against 
the angular separation and the contrast ratio axes. The exo¬ 
planet’s orbit is therefore displayed in a separation-contrast 
phase diagram. The wavelength is the V-band. The “x” indi¬ 
cates the position of the exoplanet at an arbitrary observation 
time. The ‘diamond indicates the position of the planet one 
sixth of the orbital period '/’later (giving the forward direction 
of time). The 10 plus signs indicate the evolution of the planet 
along the curve; they are linearly spaced in time. The empty and 
filled triangles are periapsis and apoapsis: r mm = 0.026,4 U and 
r max = 0.0544 6 '. The empty and filled circles are minimum and 
maximum phase angle: or m ,„ =75° and a max = 105°. The bot¬ 
tom left panel shows the orbit as seen by the observer, where the 
shaded area corresponds to the portion of the orbit that is behind 
the plane of the paper. The dashed lined corresponds to the limit 
of the same instrument sensitivities as in the preceding example. 
The right-hand side y-axis is given in photons per VLT collect¬ 
ing area per hour integrated over V-band. Two successive plus 
signs correspond to 7h46min. 



It shows that this instrument could observe Alpha Centaury Bb 
during » 35% of its orbit if /=15°. Longer integrations could 
increase the sensitivity at the expense of larger apparent orbital 
motion on the detector, which will significantly limit the angular 
differential imaging (ADI) processing capability. 

Figure [5] shows SPOC diagramsfor the two sets of orbital 
parameters Cen a 0.34 and Cen a 0 , and for all four Z. 

This figure shows that for higher orbital inclinations, the ra¬ 
dius inferred from the orbital velocity profile decreases, which 
leads to a smaller reflecting surface and consequently, a lower 
contrast ratio. The effect of eccentricity is clearly seen: curves 
for Cen a fl all run a “back and forth locus” because of their or¬ 
bital symmetry. Eccentric orbits display an increasingly complex 
shape as the inclination increases: this is due to the projection of 
the orbit locus on the sky as an observer would see it. 

We refer to Appendix [A] for more detailed studies of diverse 
SPOC diagramsthat are based on known exoplanets. 

3.2. Separations and contrast ratios for the solar system 

An interesting illustration of the usefulness of the SPOC dia- 
gramsis provided by our solar system as viewed by an external 
observer. 
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Planet 

Diameter to Earth 

a[AU] 

e 


w p [°] 

Ag 

Mercury 

0.382 

0.3870 

0.205 

7.00 

77.45 

0.142 

Venus 

0.948 

0.7233 

0.006 

3.39 

131.53 

0.67 

Earth 

1 

1.000 

0.016 

0.00 

102.94 

0.367 

Mars 

0.532 

1.523 

0.093 

1.85 

336.04 

0.170 

Jupiter 

11.2 

5.203 

0.048 

1.30 

14.75 

0.52 

Saturn 

9.44 

9.537 

0.054 

2.48 

92.43 

0.47 

Uranus 

4.00 

19.19 

0.047 

0.76 

170.96 

0.51 

Neptune 

3.88 

30.06 

0.008 

1.76 

44.97 

0.41 


Table 5. Planetary and orbital parameters for the planets of the solar system, in the V-Band. Source: http: //nssdc. gsfc. nasa. 
gov/planetary/factsheet/ Inclinations is this table are measured from the plane of the ecliptic. 



1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 


Fig. 4. Same caption as Figure [ 3 ] except that it represents the 
apparent orbital motion in mas per hour as gray scale, plotted 
against the angular separation and the polarized contrast ratio 
axes. The red dotted line corresponds to the SPHERE-ZIMPOL 
performance for 1 hour exposure (observability domain being 
above the line). 



Fig. 5. Contrast ratio versus angular separation plots for both or¬ 
bital parameters of Alpha Centaury Bb as described in Figure [3] 
and for inclinations in [15, 35, 60, 85]°. Dotted lines correspond 
to e-0 and « p =0°and full lines to e=0.34 and u> p =246°. A color 
version is available in the online journal. 


Figure [6] gives the contrast-separation diagram for the plan¬ 
ets of the solar system, whose parameters are listed in Table [5] 
In this figure, polarization is not shown (for clarity). It displays 
values for two different angles of inclination for the observer, 
0°and 45°. The observational wavelength is visible light. Note 
that no rings were modeled for Saturn; they are anticipated to 
have a variable effect on both contrast ratio and polarization (see 
discussion in Sect. 3.3). 

Note that the angular separation of Earth to the Sun is 
~ larcsec when the ecliptic inclination is 0°, consistent with 
the definition of a parsec. Earth’s average angular separation de¬ 
creases as the orbit appears more inclined to the observer. 

Figure[6]highlights the strong dependency of the contrast ra¬ 
tio and angular separation on inclination, argument at periapsis, 
and eccentricity. Mercury especially, but also Mars and Saturn 
do not follow a “back and forth locus” because the symmetry of 
their orbit as it appears to the observer is broken. 

Venus is seen to be the second brightest planet in the solar 
system (for a ringless Saturn). Significantly brighter than Earth, 
but fainter than Jupiter. A rule of thumb can explain this eas¬ 
ily: Venus receives twice as much solar flux as Earth and has an 
albedo nearly twice as high. Because the diameters of the two 
planets are close to each other, they present almost the same re¬ 
flective area (Venus has =90% of that of Earth): Venus is there¬ 
fore on average nearly four times brighter than Earth to an ex¬ 
ternal observer. The same rules can be applied to compare the 
relative brightness of Jupiter and Venus’: Jupiter is only twice 
as bright as Venus to the observer. The same rule of thumbs can 
explain the surprising faintness of Mars, Uranus, and Neptune. 
They show a relatively small reflective area because of their dis¬ 
tance to the Sun, which makes them receive little flux; further¬ 
more, Mars has a significantly low albedo. Note that these com¬ 
parative values are averaged over the entire orbits. 

Finally, this same figure highlights once more that focusing 
the discussion on averaged contrast ratios and separation is fu¬ 
tile: depending on the respective configuration of the planets. 
Earth might easily become the brightest planet in the solar sys¬ 
tem at some epochs. As an example, let us consider the case 
when the observer sees the Solar System with an inclination of 
45°. Figure [7] shows the contrast ratios of Jupiter, Venus, and 
Earth as a function of time during a whole Jovian orbit. 

It can be calculated that, during the Jovian revolution (11.9 
years): 

- CR(Earth) 6 CR(Venus) for 2.55 years (21.53% of T Jup ), 

- CR(Earth) 6 CR(Jupiter) for 1.04 years (8.82%), 

- CR(Venus) 1 CR(Jupiter) for 3.16 years (26.61%), 

- Venus is the brightest of the 3 planets for 3.11 years 

(26.16%), 

- Earth is the brightest of the 3 planets for 4.83 months 
(3.39%), 
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Angular Separation [mas] 


Fig. 6. Contrast-separation diagram for the planets of the solar system as seen from lpc distance in visible light. Small black curves 
correspond to an inclination of the ecliptic of 0°(observer reference frame); colored curves show the solar system at an inclination 
of 45°. Note that both axes have logarithmic scales. A color version is available online. 



Fig. 7. Evolution of the contrast ratio to the Sun for Venus, 
Jupiter and the Earth, during an entire Jovian revolution (11.9 
years). 


- Jupiter is the brightest of the 3 planets for 8.36 years 
(70.44%), 

where CR stands for contrast ratio. These values become 
28.28%, 18.89%, 32.14%, 30.30%, 4.83%, and 64.86% respec¬ 
tively when the observer sees the solar system at an inclination of 
75°. Note also that for an inclination of 45°, Mercury is brighter 


than Jupiter during 1.75% of Tj up . This duration increases to 
11.62% for an inclination of 15°. 

Despite the large variability in apparent brightness, it re¬ 
mains true that all solar system planets are extremely challeng¬ 
ing: Earth’s relative magnitude to the Sun is about * +23mag, 
at best. Gas giant planets at the snow line such as Jupiter do not 
always offer dramatic gains in observability, as might naively be 
thought. 


3.3. Effect of extensive planetary rings 


Rings were not included in this model. To do so would at least 
double the number of parameters describing a planet (exoplanet 
obliquity, obliquity at periapsis, ring radii (inner and outer), re¬ 
flectance parameters, etc). The ring orientation to the observer 
is critical for determining reflected light from the host-star. 
Beyond the model of our own solar system, there is little ob¬ 
servational data to constrain the speculative range of ring prop¬ 
erties. However, their presence around an exoplanet would sig¬ 
nificantly change the observable properties of planets. Rings can 
either act as reflectors with a potential polarized enhancement 
or obstruct of the exoplanet illumination or irradiance, depend¬ 
ing on their apparent inclination to the star, the observer and 
their polarized-reflectance characteristics. Both reflection and 
obstruction effects are coupled for extensive systems of rings 
that project shadows onto the exoplanet atmosphere. We refer 
to Dyudina et al. (2005 1 or |Arnold & Schneider (2006 1 for a 
more detailed discussion of the impacts of rings on light curves. 
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3.4. Distinguishing radius and albedo with polarization 


Knowing both the contrast ratio (or absolute flux) of the planet 
and its polarization degree, preferably at several orbital posi¬ 
tions, allows the observer to distinguish between several atmo¬ 
spheric models of the atmosphere and finally determine its ra¬ 
dius and albedo. For a given exoplanet atmosphere and the same 
unpolarized flux measured, a planet with large radius and low 
albedo will exhibit a higher polarization degree than a small 
exoplanet with high albedo. We refer to the relation between 
peak polarization and single-scattering albedo in Madhusudhan 
|& Bu rrows (2012 >. Such radius-albedo measurements were car¬ 
ried out in Herdyiigina et ah]|2008| and [Berdyugina et aL| ( [201 lj l. 


4. Conclusions 

4.1. Model and calculation code 

We have constructed a model to fully describe the emerging ra¬ 
diation field from an exoplanet with given orbital geometry and 
reflectance parameters. This model aclculates the three main ob¬ 
servables that are relevant for direct imaging of an exoplanet: 
polarization, contrast ratio, and angular separation as a function 
of date. A key strength of this model is the relatively restricted 
number of free input parameters despite the complex processes 
addressed. Only 3 of these input parameters - eccentricity, incli¬ 
nation, and argument of periapsis - are mandatory for perform¬ 
ing a first assessment of the variability over time of the exoplanet 
observability. 

The calculation code developed is available for pub¬ 
lic use and collaborative improvement on the python pack¬ 
age index https://pypi.python.org/pypi under “ex- 
ospoc”, together with its documentation. This code is imple¬ 
mented on the http://exoplanet.eu/ exoplanet database 
(see Schneider et al. ( |2011] >), where SPOC diagramsare accessi¬ 
ble in the exoplanet information sheets at the link “Observability 
Predictor”. 

4.2. SPOC diagrams 

We have introduced a novel tool - the SPOC diagram. It high¬ 
lights the interdependency of polarization, contrast ratio, and an¬ 
gular separation to the host star for an input exoplanet and gives 
the observer critical information for predicting expected exo¬ 
planetary signal from a minimum set of parameters. The com¬ 
plex shape of the SPOC diagramcurves highlights the fact that 
in the general case (inclination higher than :«10°), the critical 
observables strongly depend on the geometry of the orbit as it 
appears to the observer. As a consequence, we stress that re¬ 
stricting consideration to values averaged along the entire orbit 
for angular separation, contrast ratio, or polarization may be a 
misleading oversimplification in many applications. An illustra¬ 
tion of this is that an external observer would report Earth to 
be the brightest of the eight planets for a significant fraction of 
random observations. 

The variation of the main observational criteria is critical 
in timing a direct observation of an exoplanet. This variation 
mainly relies on the combined effect of the inclination, the ec¬ 
centricity and the argument of periapsis, hence the shape of the 
orbit as it appears to the observer. This latter parameter has a 
surprisingly important role to play in the visibility prediction for 
a planet in acting to either amplify or cancel the effects of incli¬ 
nation and eccentricity over the contrast ratio and angular sep¬ 
aration values over time. The albedo parameters only shifts the 


lowest and highest values of the contrast ratio; however they do 
not significantly change the span of its minimum and maximum. 
The planet radius, semi-major axis and, observer-host-star dis¬ 
tance parameters simply scale all observables to higher or lower 
contrast ratios and angular separations. 

The simultaneity of the maxima of three main observational 
criteria can also be studied with SPOC diagrams. The simultane¬ 
ity of contrast ratio and polarization (or of contrast ratio and an¬ 
gular separation) maxima is unlikely (or very unlikely) in the 
tcR.max ± (20% of T) temporal window. Therefore, the best-case 
scenario for detection, maxima of contrast, separation, and po¬ 
larization, is extremely unlikely. 
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Appendix A: Observability predictor examples 



Fig.A.l. Contrast-separation diagram for V-band reflectance pa¬ 
rameters (no polarization information added for clarity reasons) 
showing the famous HR 8799 exoplanetary system with its four 
confirmed candidates. The planet e was assumed to be the same 
size as the three other planets (1.2 Rjupiter)- All four planets but 
planet d follow a “back and forth” locus because of their null ec¬ 
centricity. Despite their large reflecting area, the contrast ratios 
are extremely low because of their distance to the host star. 



Fig. A.2. SPOC diagram for HD 80606 b, for which all orbital 
parameters are known down to a very high accuracy. The very 
high eccentricity (0.93) of this transiting planet gives a peculiar 
phase diagram. The higest contrast ratio is not reached at periap- 
sis because of the combination of i, e, and orbital parameters. 
Photon count is given per hour per VLT collecting area. The re¬ 
flectance parameters were assumed to be Jupiter-like, in V-band. 
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